Formation of naked singularities in five-dimensional space-time 
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Abstract 

We numerically investigate the gravitational collapse of collisionless particles in spheroidal configurations both in four and 
five-dimensional (5D) space-time. We repeat the simulation performed by Shapiro and Teukolsky (1991) that announced an 
appearance of a naked singularity, and also find that the similar results in 5D version. That is, in a collapse of a highly prolate 
spindle, the Kretschmann invariant blows up outside the matter and no apparent horizon forms. We also find that the collapses 
in 5D proceed rapidly than in 4D, and the critical prolateness for appearance of apparent horizon in 5D is loosened compared 
to 4D cases. We also show how collapses differ with spatial symmetries comparing 5D evolutions in single-axisymmetry, SO(3), 
and those in double-axisymmetry, U(l)xU(l). 

PACS numbers: 04.20.Dw, 04.20.Ex, 04.25.dc, 04.50.Gh 



I. INTRODUCTION 

The so-called "large extra-dimensional models" as a 
consequence of brane-world pictures have changed our 
viewpoints for a way of understanding the fundamen- 
tal forces. The scenarios of unifying the gravity at 
TeV scale or so open the possibility of verification 
of higher-dimensional space-time models at the CERN 
Large Hadron Collider (LHC). If the LHC detects pro- 
ductions (and evaporations) of mini black-holes as ex- 
pected, then humankind will encounter a Copernican 
change of our outlook of the Universe. 

With this background, black-holes in higher dimen- 
sional space-time are extensively studied for a decade. 
Many interesting discoveries of new solutions have been 
reported, and their properties are being revealed. How- 
ever, fully relativistic dynamical features, such as the for- 
mation processes, stabilities and late-time fate, are still 
unknown and they are waiting to be studied. Several 
groups including us begin reporting numerical studies in 
various topics in higher-dimensional models [H-Q- 

We, in this article, report our numerical simulations 
on gravitational collapse in axisymmetric space-time. 
The topic has been studied in many ways in (3 + 1)- 
dimensional space-time (4D, hereafter), among them the 
most impressive result we think is the work by Shapiro 
and Teukolsky [?J (ST91, hereafter); a highly prolate 
matter collapse which may form a naked singularity. We 
repeat their simulations and also compare them with 
(4 + l)-dimcnsional (5D) versions. 

In classical general relativity, it is well known that a 
space-time singularity will be generally formed in gravita- 
tional collapse of non-singular asymptotically flat initial 
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data. If a singularity forms without an event horizon, all 
physical predictions become invalid. In order not to occur 
such a disastrous situation, Penrose proposed the cosmic 
censorship conjecture [8|, which states that singularities 
are always clothed by event horizons. 

On the other hand, for nonspherical gravitational col- 
lapses, Thorne proposed the hoop conjecture [9j which 
states that black holes with horizons are formed when 
and only when a mass gets compacted into a small re- 
gion. He expressed the compactness with a 'hoop ' around 
matter. If matter configuration is highly aspherical, then 
the hoop length becomes larger. If so, the conjectured 
inequality does not hold, i.e. horizon will not be formed. 
Thereby, the hoop conjecture indicates that a highly as- 
pherical matter collapse will lead to a naked singularity. 

ST91 numerically showed that axisymmetric space- 
time with collisionless matter particles in spheroidal dis- 
tribution will collapse to singularity, and there are no ap- 
parent horizon formed when the spheroids are highly pro- 
late. The behaviors are consistent with their initial data 
analysis and support the hoop conjecture. How- 

ever, since numerical evolutions cannot provide the final 
structures nor the conclusive information for formation 
of naked singularities, debates were raised after their an- 
nouncement. For example, Wald et al [ll], UH showed 
examples that 3-dimensional hypersurfaces can hit sin- 
gularities nevertheless the space-time is consistent with 
cosmic censor. Their examples are not directly related 
with the numerical results with ST91, but we learned 
that a numerical result provides only a limited evidence. 

Regarding to the 5D cases, the hoop conjecture is sup - 
posed to be replaced with the hyper-hoop version jl3l4l7| . 
i.e. a criteria is not a hoop but a surface. In our previ- 
ous work[IJ, we numerically constructed initial data se- 
quences of non-rotating matter for 5D evolutions and ex- 
amined the hyper-hoop conjecture using minimum area 
around the matter. We found that the areal criteria 
matches with the appearance of apparent horizons for 
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spindle matter configurations (but not for ring configu- 
rations). The sequences suggest that a highly prolate 
spindle in 5D will form a naked singularity similar to the 
4D cases. We also found that a condition for a naked 
singularity formation is relaxed in 5D than 4D cases. 

One of the objectives of the present work is to com- 
pare the dynamics between 4D and 5D. A simple estima- 
tion from the free-falling time indicates that the gravi- 
tational collapse in 5D takes longer time than 4D cases. 
We show this is not applicable to the highly non-linear 
final stages. In 5D, two axes can be settled as rota- 
tional symmetric axes, so that we also compare gravita- 
tional collapses in axisymmetry with those in "doubly" - 
axisymmetric space-time. 



II. NUMERICAL CODE 

We evolve five-dimensional axisymmetric [symmet- 
ric on z-axis, SO(3)] or doubly-axisymmetric [symmet- 
ric both on x and z-axes, U(l) xU(l)], asymptotically 
flat space-time (see Figure [lj. For the comparison, 
we also performed four-dimensional axisymmetric space- 
time evolutions. 

We start our simulation from time symmetric and con- 
formally flat initial data, which are obtained by solving 
the Hamiltonian constraint equations |jj . The asymptot- 
ical flatness is imposed throughout the evolution, which 
settles the fall-off condition to the metric as ~ 1/r for 
4D cases and ~ 1/r 2 for 5D cases. 

The matter is described with 5000 collisionless par- 
ticles, which move along the geodesic equations. We 
smooth out the matter by expressing each particle with 
Gaussian density distribution function with its typical 
width is twice as much as the numerical grid. The parti- 
cles are homogeneously distributed in a spheroidal shape, 
parametrized with a and b (Figure [J), or eccentricity 
e=^/l-a 2 /b 2 . 

By imposing axisymmetry or double-axisymmetry, our 
model becomes practically a (2 + l)-dimensional problem. 
We construct our numerical grids with the Cartesian co- 
ordinate (x, z), and apply the so-called Cartoon method 
2, 18] to recover the symmetry of space-time. 
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FIG. 1: We evolve five-dimensional (a) axisymmetric [SO(3)] 
or (b) double-axisymmetric [U(l)x U(l)], asymptotically flat 
space-time using the Cartesian grid. The initial matter con- 
figuration is expressed with parameters a and b. 



The space-time is evolved using the Arnowitt-Deser- 
Misner (ADM) evolution equations. It is known that 
the ADM evolution equations excite an unstable mode 
(constraint- violation mode) in long-term simulations (l9l . 
|20|. However, we are free from this problem since gravi- 
tational collapse occurs within quite short time. By mon- 
itoring the violation of constraint equations during evo- 
lutions, we confirm that our numerical code has second- 
order convergence, and also that the simulation contin- 
ues in stable manner. The results shown in this report 
are obtained with numerical grids, 129x129x2x2. We 
confirmed that higher resolution runs do not change the 
physical results. 

We use the maximal slicing condition for the lapse 
function a, and the minimal strain condition for the shift 
vectors (3 l . Both conditions are proposed for avoiding the 
singularity in numerical evolutions 21J, and the behav- 
ior of a and /3 l roughly indicates the strength of gravity, 
conversely. The iterative Crank-Nicolson method is used 
for integrating ADM evolution equations, and the Runge- 
Kutta method is used for matter evolution equations. 

For discussing physics, we search the location of appar- 
ent horizon (AH), calculate the Kretschmann invariant 
(I = RabcdR ahcd ) on the spacial hypersurface. 



III. RESULTS 

We prepare several initial data fixing the total ADM 
mass and the eccentricity of distribution, e = 0.9. By 
changing the initial matter distribution sizes, we observe 
the different final structures. Figure [5] shows snapshots 
of 5D axisymmetric evolutions of model b/M = 4 and 
10 (model 5DS/3 and 5DS5, respectively; see Table HJ; the 
former collapses to a black hole while the latter collapses 
without AH formation. 

All the models we tried result in forming a singularity 
(i.e., diverging I). We stopped our numerical evolutions 
when the shift vector was not obtained with sufficient 
accuracy due to the large curvature. For model 5DS<5, 
we integrated up to the coordinate time t/M = 15.4 and 
the maximum of the Kretschmann invariant I max became 
0(1000) on z-axis (see Figure|3]), but AH was not formed. 

When the initial matter is highly prolate, AH is not 
observed. This is consistent with 4D cases 0, E3] , and 
matches with the predictions from initial data analysis in 
5D cases 0, H3 |. The location of I ma x is on z-axis, and 
just outside of the matter [22^. This is again the same 
with 4D cases Q- The absence of AH with diverging I 
suggests a formation of naked singularity in 5D. 

In order to compare the results with 4D and 5D, we 
reproduced the results of ST91. We then find that the 
e = 0.9 initial data with b/M = 10 (model 4D<5) col- 
lapses without forming AH, and the code stops at the 
coordinate time t = 20.91 with I max — 84.3 on the pr- 
axis {z/M = 6.1); all the numbers match quite well with 
ST91. (Note that our slicing conditions and coordinate 
structure is not the same with ST91.) 
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FIG. 2: Snapshots of 5D axisymmetric evolution with the ini- 
tial matter distribution of b/M — 4 [Fig.(al) and (a2); model 
5DS/3 in Table U and 10 [Fig.(bl) and (b2); model 5DS5]. We 
see the apparent horizon (AH) is formed at the coordinate 
time t/M = 3.3 for the former model and the area of AH 
increases, while AH is not observed for the latter model up 
to the time t/M = 15.4 when our code stops due to the large 
curvature. The big circle indicates the location of the max- 
imum Kretschmann invariant I max at the final time at each 
evolution. Number of particles are reduced to 1/10 for figures. 




FIG. 3: Kretschmann invariant I for model 5DS5 at t/M = 
15.4. The maximum is 0(1000), and its location is on z-axis, 
just outside of the matter. 



TABLE I: Model-names and the results of their evolutions 
whether we observed AH or not. The eccentricity e of the 
collapsed matter configurations is also shown; gah and ef are 
at the time of AH formation (if formed), and on the numeri- 
cally obtained final hypersurface, respectively. 



Figure 2] is for the comparisons of hypersurfaces for 
the 5D models which collapses (a) with forming AH, (b) 
without forming AH, and (c) 4D collapses without form- 
ing AH. We see hypersurfaces are bending due to the 
slicing conditions, and figures tell us how numerically in- 
tegrated region covers the space-time. 

We also performed 5D collapses with doubly- 
axisymmetric [U(l)xU(l)] space-time. The matter and 
space-time evolve quite similar to 5D and 4D axisymmet- 
ric cases, but we find that the critical configurations for 
forming AH is different. Table U summarizes the main 
results of 4D and two 5D cases. We find that AH in 5D 
is formed in larger b initial data than 4D cases. This re- 
sult is consistent with our prediction from the sequence 
of initial data [l| . AH criteria with initial b is loosened 
for 5D doubly- axisymmetric cases. 

We also show the eccentricity of matter, eAH and ef, 
at the time of AH formation (if formed) and at the final 
time in the simulation, respectively. The numbers in Ta- 
ble U indicate that the eccentricity itself is not a guiding 
measure for AH formation, but they give us a hint for 
understanding the differences. 

In 4D, the eccentricity increases after AH is formed 
(4Da), while it decreases in 5D axisymmetric cases (5DSq 
and 5DS/3). That is, the 5D collapses proceed towards 
spherical configurations. This fact would be explained 
by the degree of freedom of the movements. In general, 
5D spacetime is expected to produce much gravitational 
radiation than 4D spacetime 0, [E[ , since there are more 
modes of oscillation exist. Gravitational radiation nor- 
mally works to change shapes to spherical, because it 
is produced from the acceleration of spacetime and car- 
ries the energy away. (It is known that compact binary 
system will evolve into a circular orbit due to the emis- 
sion of gravitational radiation.) Therefore collapses in 
5D spacetime are likely to evolve towards more spheri- 
cal. This interpretation together with the hoop conjec- 
ture will explain why AH-formation condition is loosened 
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FIG. 4: The snapshots of the hypersurfaces on the z-axis in the proper-time versus coordinate diagram; (a) model 5DS/3, (b) 
model 5DS<5, and (c) model 4D<5. The upper most hypersurface is the final data in numerical evolution. We also mark the matter 
surface and the location of AH if exist. The ranges with X > 10 are marked with bold lines and peak value of T express by 
asterisks. 



in 5D cases. 

In the 5D doubly-axisymmetric cases, on the other 
hand, the magnitude of ef is smaller than 5D axisym- 
metric cases for small b/M cases (5DSa vs. 5DUa, or 
5DS/3 vs. 5DU/3), while it is larger for large b/M cases 
(5DS7 vs. 5DU7, or 5DS<5 vs. 5DU(5). We think this is 
because that the doubly-axisymmetric collapses proceed 
in more symmetric manner than axisymmetric collapse 
near the origin, while they proceed in more 4D-like ax- 
isymmetric collapses near the axes far from the origin. 
The collapses of small b/M initial data, therefore, will 
evolve into more spherical shape, while the large b/M 
initial data will evolve increasing the eccentricity, where 
the latter is similar to 4D cases. 
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FIG. 5: The Kretschmann invariant I at the location of I max 
on the final hypersurface is plotted as a function of proper 
time at its location. Labels indicate model-names in Table U 
The time of AH formation (t=0.6 for model 5DS/3, t=0.9 for 
5DU/3) is shown by a dot. 

In Figure [H we plot X at the point which gives X m ax on 
the final hypersurface as a function of proper time. The X 
diverges at the end of simulations in all the cases, but the 
diverging time becomes later for larger b/M initial data. 
We see that 5D-collapses are generally proceeding more 
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rapidly than 4D collapses. We also see that collapses in 
5D doubly-axisymmetric space-time is proceeding more 
slowly than 5D single axisymmetric cases. If we observe 
further, the model 5DU/3 evolves quite similar to 5DS/3, 
while 5WS evolves quite similar to 4D<5. These behaviors 
support the previous understandings of the evolution of 
the eccentricity. 



IV. DISCUSSIONS 

In this article, we reported our numerical study of 
gravitational collapses in 5D space-time. We collapsed 
spheroidal matter expressing with collisionless particles, 
and observed the evolution of the Kretschmann invariant 
and the apparent horizon (AH) formation. 

The collapsing behaviors are generally quite similar 
to the cases in 4D, but we also found that (a) 5D- 
collapses proceed rapidly than 4D-collapses, (b) AH ap- 
pears in highly prolate matter configurations than 4D 
cases, (c) doubly-axisymmetric [U(l)xU(l)] assumption 
makes collapse proceed towards more spherical when it 
forms AH, but presents quite similar behavior with 4D 
cases for large configurations, and (d) the positive evi- 
dence for appearance of a naked singularity in 5D. 

Up to this moment, we only checked the existence of 
apparent horizons, and not the event horizons. The sys- 
tem does not include any angular momentum. We are 
implementing our code to cover these studies. 

We are now preparing our next detail report including 
the validity of hyper-hoop conjecture in 5D, and the cases 
of the ring objects. 
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